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Abstract. In the present short communication, we revisit nonlinearity management of the time- 
periodic nonlinear Schrodinger equation and the related averaging procedure. We prove that the 
averaged nonlinear Schrodinger equation does not support the blow-up of solutions in higher 
dimensions, independently of the strength in the nonlinearity coefficient variance. This conclusion 
agrees with earlier works in the case of strong nonlinearity management but contradicts those 
in the case of weak nonlinearity management. The apparent discrepancy is explained by the 
divergence of the averaging procedure in the limit of weak nonlinearity management. 
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1. Introduction 

In the past few years, there has been a large volume of literature regarding the applications of 
the nonlinear Schrodingcr equation (NLS) in the presence of the so-called nonlinearity management 
(often referred to also as Feshbach resonance management). The NLS is a prototypical dispersive 
nonlinear wave equation of the form: 

iu t = -Au + T(t)\u\ 2 u + V(x)u, (1) 

where u(x,t) is a complex envelope field, V(x) > is an external potential, T(t) is a time-periodic 
nonlinearity coefficient, and A is the Laplacian operator with x € R , d > 1. 

Nonlinearity management arises in applications in optics for transverse beam propagation in 
layered optical media PP, as well as in atomic physics for the Feshbach resonance of the scattering 
length of inter-atomic interactions in Bosc-Einstcin condensates (BECs) 2 . In the latter case, the 
periodic variation of T(t) through an external magnetic field has been used as a means of producing 
robust matter- wave breathers in quasi-one-dimensional BECs 0- It has also been suggested 
that the nonlinearity management may prevent collapse-type phenomena in higher dimensions 4 . 
Theoretical studies of the nonlinearity management were performed with a reduction of the time- 
periodic PDE problem to a time-periodic ODE problem using a variational method [5] and a 
moment method [2]. 

The physical relevance of the time-periodic NLS equation, as evidenced by the above works, 
led to further developments in analysis of the PDE problem JJJ. As an example, the success of 
the averaging theory for optical solitons in the presence of strong dispersion management, led 
to an analogous development for strong nonlinearity management produced originally in [S] and 
systematized in The time-periodic NLS equation JJJ is averaged in the limit e — * 0, where e 
measures a short period and the large variation of T(t) in the decomposition 

r = 70 + \l (-) , (2) 

where j(t), t = t/e has a unit period and zero mean. After the averaging procedure, the time- 
periodic PDE problem is governed by the averaged NLS equation, 

iw t = -Aw + j \w\ 2 w + V(x)w - a 1 (|V|w| 2 | 2 + 2|w| 2 A|w| 2 ) w, (3) 

where ^ 

a 2 = [ -y^Mdr, (4) 
Jo 

and 7_i(t) is the mean-zero anti-derivative of 7(t). Derivation and local well-posedness of solutions 
of the averaged NLS equation © in function space H°°(M.) are reviewed for d = 1 in j^j. 

In the present paper, we consider global well-posedness of solutions of the averaged NLS 
equation @ in the energy space i? 1 (R d ) for d > 1. The use of H 1 and d > 1 seems more appropriate 
for physical applications of the averaged model JSJ- In particular, we address the question whether 
the averaged NLS equation with a nonlinearity management (a 2 > 0) arrests the blowup of solutions 
of the NLS equation in two and three dimensions that would occur if no nonlinearity management 
was applied (er 2 = 0). We show that the averaged NLS equation is globally well-posed and no 
blowup of solutions occurs for a 2 > 0. This is demonstrated initially, in Section 2, from the point of 
view of scaling arguments. The statement is subsequently proved using rigorous estimates in Section 
3. In Section 4 we compare the above conclusion and earlier works where possibilities of blowup of 
solutions of the full time-periodic NLS equation Q have been reported. Section 5 summarizes our 
findings. 
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2. Formal Scaling Arguments 

The averaged NLS equation has a standard Hamiltonian form (see [0]) with the Hamiltonian 
functional: 

H(w) = J (|Vw| 2 + !^M 4 + V{x)\w\ 2 + a 2 \w\ 2 |V|w| 2 | 2 ) dx. (5) 

R d 

Due to the gauge invariance, the averaged NLS equation @ also conserves the squared L 2 norm: 

P(w) = I \w\ 2 dx. (6) 



Solitary wave solutions of the averaged NLS equation © are critical points of H(w) at the level set 
of fixed values of P(w). 

Using formal scaling arguments JUj (see also the review of Ref. JI]), we consider a two- 
parameter family of dilatations: 

w = bW(ax), (7) 

where (a, b) are parameters and W(£) is a suitable function of £ = ax. The squared L 2 norm JBJ is 
preserved by the dilatations whenever b — a d l 2 . The Hamiltonian (JSJ at the dilatations (7J is 
scaled as a function of parameter a > 0: 

ff(o) = Jo(a) + a 2 Ji + To^Ia + <J 2 a 2d+2 I 3 , (8) 

where 



/i = j \^W\ 2 dt h = \J \W\ A d^ h = J \W\ 2 \V\W{< 



and 

Io(a) 



J V (?) \W\ 2 d£. 



Let us consider the case of no nonlinearity management and no external potential, when a 2 = 
and V(x) = 0. It follows from © that the Hamiltonian function H(a) is positive definite in the 
dcfocusing case, when 70 > 0. In the focusing case, when 70 < 0, the Hamiltonian function H(a) 
is bounded from below for d = 1 and d — 2, 7 cr < 70 < and is unbounded from below for d = 2, 
70 < 7 cr and d = 3, where 

h 
H 

When H(a) is unbounded from below as a — > 00, the critical points of H(w) at a fixed value of 
P(w) (i.e., solitary wave solutions) can not be stable for small width a -1 and instability of solitary 
waves implies a blowup of localized initial data in the time evolution of the cubic NLS equation (see 
[III for details). 

When the nonlinearity management is applied, the last term in the decomposition (jHJ always 
dominates and it preserves the boundness of H(a) from below for any a 2 > 0. This indicates on 
the level of formal scaling arguments that the blowup of solutions is arrested by the nonlinearity 
management term in the averaged NLS equation We shall prove this conjecture with rigorous 
analysis of well-posedness of solutions. We also note that the first term in the decomposition (jHJ does 
not change the conclusions above if V(x) is a smooth non-negative potential, such that Io(a) > 0. 
Typical examples of V(x) are parabolic magnetic traps, when V ~ x 2 , and periodic optical lattices, 
when V ~ sin 2 (fco^)- 
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3. Rigorous Analysis of Well-Posedness 

The rigorous analysis of the local well-posedness of the averaged NLS equation is not a trivial 
task. In fact, to the best of our knowledge, one cannot verify even the local existence and uniqueness 
of solutions to this problem. The problem has been considered in one dimension d = 1 by Poppenberg 
where a local well-posedness result with data in H°° was established |12| . In higher dimensions 
d > 2, one needs to require that the initial data be in H s ' m for sufficiently large s,m, that is 
J R „(l + |a;|) 2m (|w;o| 2 + |(9 s wo|) 2 d2; < oo (see ^3]). In addition, one needs to assume a "non-trapping" 
condition on the symbol of the second order operator. This is a geometric condition, which depends 
on the profile of the initial data (see for details). 

To summarize, we cannot state a precise condition under which the averaged NLS equation © 
has a (local) solution that preserves values of P and H constant in time. We will however show, 
that whenever such a local solution exists for a short time < t < to < oo, it can be extended 
globally for all t > 0. 

To that end, we represent the Hamiltonian H(w) in the form: 



where 



and 



H(w) = Hi(w) +j H 2 (w), 
Hx(w) = J (\Vw\ 2 + V(x)\w\ 2 + a 2 \w\ 2 \V\w\ 2 \ 2 ) dx > 



1 ' ■ ,4, 



H 2 (w) = -J \w\*dx > 0. 

R d 

We consider the focusing case 70 < and prove that Hi(w) and Pi2(w) are bounded by the two 
conserved quantities H(w) and P(w). 

First, we quote a variant of the Gagliardo-Nirenberg inequality (see (1.1.16) on p. 15 in |14p. 

Lemma 1. For all 1 < p,q,r < 00, 9 £ (0, 1) and r _1 = Op^ 1 + (1 — 9)q~ 1 , it is true for every 
function f(x) on x £ M. d that 

ll/ll^ < ll/lD/lli^ (9) 

Next, we define and use the Fourier transform and its inverse for a function f{x) on x € R d : 
/(£) - / f{x)e- 2 "^dx, f(x) = / /(Oe^'flde. (10) 



The Plancherel's formula gives ||/|| L 2 = ||/||l 2 and the inequality ||/|| L oo < H/Hl 1 is immediate 
from the definition. 

Let x( s ) be a C X) (R+) function with a compact support on s G [0,2], such that x( s ) = 1 
on s £ [0,1]. Let Xd( x ) be a Co°(R d ) function, such that Xd( x ) — x(\ x \)- The (smooth) Fourier 
multiplier P<jy is defined for every positive N by P<Nf{£.) — Xd(£,/N) /(£). Equivalcntly, 

P <N f(x) =N d [ Xd(N(x - y))f{y)dy. (11) 



Let PyN = Id — P<at. Then both P>n,P<n are self-adjoint, bounded on L 2 operators. The 
following statement is a modification of the Sobolev embedding theorem. 

Lemma 2. There exists a constant Cd > 0, which depends only on the dimension d > 1, so that it 
is true for every function fix) on x £ M. d that 

\\f\\ L *<C d i\\Vf\\ L2 + \\f\\ L1 ). (12) 
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Proof. Since / = P<if + P>if, it will suffice to show that 

\\P<if\\ L *<C d \\f\\ Li (13) 
||iW|| L 2 <C||V/|| L2 (14) 

The bound (|14|l follows simply by the Plancherel's theorem and — 27ri£j/(£). Indeed, 



p>i/iir.< / i/(0i 2 ^< / iei 2 i/(0i 2 ^ 



4tt 2 



|V/(0| 2 dC < 



k c l>i 



|€|>i 



l«l>i 



<^ IWmM = — 2 \\Vf\\ 2 L2 



4tt 2 ' 



By duality, the bound (| 1 31) follows from the estimate ||-P<i/||^oo < Cdll/llxa- That is trivial as well, 
since 



\P<if\\ L °° < 



P<if 



< 



ki<2 



\mm<c d ( j \mM) l/2 <c d \\f\\ L2 . 

I«l<2 



The positive constant C d can be taken to be the square root of the volume of the ball in R d with 
radius 2. □ 

The central result of our analysis is the following theorem. 

Theorem 1. There exist e(d) > and C(e, d, a) > 0, so that for every < e < s(d), it is true for 
every function 4>{x) on x £ M. d that 



i4 

L 4 



<eH 1 (^)+C(e,d J <r)(w 
Proof. Let / = h 3 / 2 with h(x) > in (JT^J and obtain 

1/2 



2 

L 2 



ft d da; < Ci 



Next, we set h = \cf>\ 2 and obtain 




h\Vh\ 2 dx 



10111=) 



/i 3/2 cte 



(15) 



| Wdx<C 2 | || H 2 |V|0| 2 | 2 ^ 

2 I |0| 2 |V|0| 2 | 2 ^ + 2 || |^| 3 dx 




|</>| 3 dx 



< c d , a (Hi(<i>) + H\\ls), 

for some positive constant Cd )(T . We have used here that 

(Va~ + b) 2 <2(a + b 2 ). 

By the Gagliardo-Nirenberg inequality JjJJ, we have ||</>|| L 3 
inequality is rewritten in the form: 



< 



\ 1/2 \\6\\ 1/2 



such that the last 



HWle < C d ,„H x {cj>) + C 2 ||0|| 3 6 ||0|| 3 2 < C^H^) + C 2 ( eUWle + -||0|| 6 L2 ) , (16) 
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where in the last line we have used the Cauchy-Schwartz inequality: 

b 2 

Ve>0: ab<ea 2 + — . (17) 
4e 

Let e < s(d), where 2C 2 e(d) = 1. Then, the term H^H^e can be estimated from the bound l|16|) as 
follows: 

Hf L 6 < c^micj)) + c d uf L2 (is) 

for some constants Cd, a > and Cd > 0. By the Gagliardo-Nirenberg inequality we have 



— H^Hl^H^Hi^ 4 ! sucn the upper bound for \\4>\\ L 4 follows from l(T7|l and l{T5)l: 



which is the desired upper bound Q15p. □ 

As a corollary of the main theorem, we pick e = e(d)/2 and immediately obtain the following 
upper bounds. 

Corollary 1. There exists constants C\ > and C*2 > i/iai depend on d,"/o,a, so that 

Hi(w) <C 1 (H{w)+P(w)+P 2 (w)) (19) 

i? 2 (w) < C 2 (H(w) + P(w) + P 2 {w)) . (20) 

Since H(w) and P(w) are conserved in the time evolution, the Cauchy problem for the averaged 
NLS equation © has global solutions in the energy space H 1 (M. d ), if the initial data w(x,0) gives 
rise to local solutions in the same energy space. Therefore, the blowup of solutions of the NLS 
equation with a 2 = in d > 2 is arrested by the nonlinearity management for any a 2 > 0. 

4. Averaged Equation versus Full Dynamics 

We have proven that the blowup does not occur in the averaged NLS equation J3J) in higher 
dimensions d > 2 for a 2 > 0. This result raises the question whether the blowup of solutions 
is arrested in the full NLS equation Q for any non-zero variance of the time-periodic nonlinearity 
coefficient L(i). We address this question within the ODE reduction of the time-periodic problem, 
which was considered recently with a variational method 4, 5 and a moment method In both 
cases, the time evolution of the radially symmetric localized solutions of the full NLS equaiton 
is approximated by a time-dependent, generalized Ermakov-Pinney |15j equation: 

m = %+nt)^, (2i) 

where R(t) > is an effective width of a localized solution, while (Qi, Q2) are constants found from 
an initial data, such that Q2 > (see 6 for details). We shall consider the critical case d = 2 and 
rewrite the ODE (|21() with the nonlinearity coefficient T(t) in in the explicit form: 

R(t) = ^±^M, (22) 

where a = Q± + 70Q2, P — Q21 'e > 0, and 7 is a mean-zero e-periodic function of t. Conditions for 
blowup and existence of bounded oscillations in solutions of the ODE ill' 21) were recently reviewed in 
jH] (see also references therein). The sufficient condition for the blowup (when R(t) — > in a finite 
time t — > to) is 

a + 8 max (7) < 0. (23) 

0<t<e 
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The sufficient condition for the unbounded growth of the solution width (when R(t) — > oo as t — > oo) 
is 

a + f3 min (7) > 0. (24) 

0<t<e 

The necessary condition for the bounded oscillations of R(t) > for any t > is 

a < 0, a + (3 max (7) > 0. (25) 

0<t<e 

Numerical simulations in [^j showed that the condition (|25(l was also sufficient in the case a < for 
the blowup arrest for any t > 0. 

It is obvious from the explicit scaling that j3 3> \a\ in the asymptotic limit e — > 0. Therefore, 
the second condition (|25|l is satisfied for a < and solutions of the time-periodic ODE (|22H do not 
collapse, in agreement with our results derived for the averaged NLS equation 

Previous work (see also the review in [5]) has also addressed the averaged NLS equation 
in the limit of weak nonlinearity management, when 7(7") is rescaled as 7 = £j{t) and the parameter 
a 2 is small in the limit e — > as a 2 = e 2 a 2 . Solutions of the averaged NLS equation J2J at the 
leading order e = blow up in a finite time, but the small e 2 -terms formally stabilize the blow-up 
for any a 2 > 0. Although this approximation of the critical NLS equation has been considered in 
many applications of nonlinear optics (see Sections 4-5 in |16| and references therein) , it is clearly 
insufficient for a correct identification of the domain, where the blowup of solutions occurs. Indeed, 
while the averaged NLS equation © with small a 2 predicts no blowup of solutions, it is clear 
that the weak nonlinearity management corresponds to the case (3 « \a\ and solutions of the ODE 
problem IL'l'l) (and those of the full PDE problem JIJ) collapse in the domain l|23|l . 

Following the work , we address the failure of the averaging procedure for weak nonlinearity 
management in a simple time-periodic ODE problem: 

m = a + PM**r) t T = i f (26) 

where a < 0, (3 > and (a, (3) are order of O(l) in the limit e — > 0. By using the formal 
asymptotic multi-scale expansion method (see 9 for details), we construct an asymptotic solution 
to the problem 126(1 : 

R = r{t)+e 2 R 2 (T 7 r)+e i R 4 (r,r) + 0(e 6 ), r=-, (27) 

where R2 and i?4 arc recursively found from the set of linear inhomogeneous problems, 

13 

3a/? . . , 3/3 2 . v 

Ri = ~ sm(27rT + 4 7 cos(47rr). 

(27r) 4 r' 8(27r) 4 r' 

The mean- value term r(t) satisfies an extended dynamical equation that excludes secular growth of 
the correction terms of the series (12711 in r: 

a , . 2 3/3 2 , . 4 15a/3 2 , ^.^ „ e 



r = -^t + e ■—=■ + e 



TT + 0(n t=-. (28) 



r 3 2r 7 2r n v n 2ir 



The averaged ODE problem (|28|l is an equation of motion for an effective particle with a coordinate 
r(t) in the potential field with an effective potential energy: 

Uir) = ^^ + e^ + { f). (29) 

When a < and e = 0, the particle with r(0) > reaches r = at a finite time t — to < 00, 
that indicates the blowup of a localized solution. When the next e 2 -term is taken into account 
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Figure 1. Numerical simulations of the full NLS equation Q (solid curve) and the averaged NLS 
equation J3J with a fourth order in time scheme, where spacings are dx = 0.075 and dt = 10 -5 . 
The half-width s of the wavefunction is shown as a function of time t. 

(as in the approximation of weak nonlinearity management 01 15]), the blow-up is arrested and the 
mean- value term r(t) oscillates in an effective minimum of the potential energy U(r), truncated 
at e 2 -terms. When the next e 6 -term is taken into account (beyond the approximation of weak 
nonlinearity management), the potential energy U(r) with a < does not prevent the blowup of 
the localized solution depending on the initial data r(0). Existence versus non-existence of blowup 
depends on the ratio of parameters (a, (3) but the difference can only be detected in the averaging 
method if convergence of the power series l|29|l is established in a closed analytical form. 

Similarly, the averaged NLS equation can not be used in the limit of weak nonlinearity 
management for an accurate prediction of existence versus non-existence of blowup of solutions. In 
order to illustrate this point, we have performed numerical simulations of the full NLS equation (JTJ 
in d = 2 with 

r(i) = -20.76 + 8sin(27rf). 

We have observed that collapse of localized initial data does occur (see solid curve on Fig. ^) by 
monitoring the half-width (of one dimensional slices along y = 0) of the wavefunction (for radially 
symmetric Gaussian initial data), 

_ 1 H x 2 \u{x,Q,t)\ 2 dx \ 1/2 
S ~ 2 ^ J \u(x,0,t)\ 2 dx ) ' 

until it becomes comparable to the lattice grid spacing used (at that scale collapse is arrested, since 
the numerical scheme cannot resolve scales below the grid spacing). On the other hand, numerical 
simulations of the averaged NLS equation J3j with the same parameters show that the half-width 
s never decreased below s < 0.23 (see dashed curve on Fig. ^) indicating the absence of collapse 
in accordance with the rigorous results presented above. Therefore, the averaged NLS equation © 
can only be used for modeling of the blowup arrest in the limit of strong nonlinearity management 
of the full NLS equation when max(7) 3> ej'yo I an d e is small. 
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5. Conclusion 

In conclusion, we have studied the global well-posedness of solutions of the averaged NLS equation 
that describe strong nonlinearity management of the time-periodic NLS equation. We have 
showed with formal scaling arguments and rigorous analysis that the blowup of solutions in 
higher dimensions is arrested within the averaged NLS equation. We have also discussed the non- 
applicability of the averaged NLS equation to the weak nonlinearity management, where the blowup 
of solutions can occur beyond the weak management limit. 

It is an open problem to study well-posedness of the full time-periodic NLS equation, depending 
on parameters of the nonlinearity management and profile of initial data. Rigorous results on the 
latter problem are only available within the ODE approximation 1|21|) . when the PDE model is 
reduced to a dynamical system with one degree of freedom. It would be particularly interesting 
to study mathematically and to examine numerically whether the theoretical prediction from the 
method of moments provides an optimal bound for the full PDE model with arbitrary initial data. 
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